A Simple Expression for the Cold Compression Curve 



V. Celebonovic 

Institute of Physics, Pregrevica 118, 11080 Zemun-Beograd, Yugoslavia 
E-mail celebonovic@exp.phy. bg. ac. yu 
and: E-mail vladanc@mrsysl.mr-net.co.yu 



Abstract: The aim of this contribution is to present expressions for the bulk modulus of a material and its 
pressure derivative obtained by using the semiclassical theory of dense matter proposed by P. Savic and R. 
Kasanin. Some possibilities for the application of these expressions are briefly discussed. 

Introduction 

The term "equation of state" (EOS) of a system denotes any kind of relationship, established experimen- 
tally or theoretically, between the parameters of state of the system under consideration. In the particular 
case of a thermomechanical system, the EOS has the general form f(p, V,T) = 0. The knowledge of this 
function is of fundamental importance in various physical and astrophysical problems (for recent examples 
see, for instance, Holzapfel, 1996; Schulte and Holzapfel, 1996; Schaab, Weber et al. 1996). 

A simplified version of the EOS, which does not take into account the influence of the temperature is 
called the cold compression curve (CCC). Strictly speaking, the CCC is a T = OK isotherm of dense matter, 
but, in practice, it can be applied at low temperatures. In theoretical studies, the form of a cold compression 
curve is a consequence of the assumptions of the underlying theory ( Holzapfel, 1996 contains an instructive 
list). One of the simplest forms of a CCC is the so-called Murnaghan EOS, which amounts to 

P(p) = (B/B')[( Po /p)- B ' -1} (la) 

In this expression, B and B' denote the isothermal bulk modulus of a material under standard conditions, 
and its first pressure derivative. These two functions are obviously material dependent, and their knowledge 
is crucial for all applications of EOS. The bulk modulus is defined as 

B = -V(dP/dV) T = P {dP/dp) T (1) 

The aim of this contribution is to propose simple expressions for the isothermal bulk modulus and its 
pressure derivative. All the calculations will be performed within a particular semi-classical theory of the 
behaviour of materials under high pressure, proposed by P. Savic and R. Kasanin in the early sixties (Savic 
and Kasanin, 1962/65). For recent reviews of their theory see Savic and Celebonovic (1994); Celebonovic 
(1995). 

The Calculations 

The necessary "ingredient" for the calculation of B and B' is the function dP/dp. In the theory proposed 
by Savic and Kasanin this function has the following form. 

dP/dp = (N A e 2 /9A)Q (2) 

where 

Q = {A/a i )f i {a i )-f' i {a i ) (2a) 
The symbol fi(cii) denotes the function 
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fi(a,i) = Ci + Bi exp(7 i z l ) (3a) 

with 

a, = (A/SNap^ 3 (3b) 

and 

Zi = (1 - a*/ ai )/(l - a; 1/3 ) (3c) 

Details of derivations of these equations are available in the original publications of Savic and Kasanin. 
The meaning of various symbols is as follows: A and Na denote the atomic mass of the material and 
Avogadro's number, a*, Bi, Ci, ji, Ui are constants (whose numerical values are known within the theory) 
in the i-th phase of a material subdued to high pressure. Physically, aj is the mean inter-particle distance 
in the i-th phase. 

Inserting eq. (3a) into (2a), it follows that 

Qi = (4/a,-)[Cj + Bi exp(7iZi)] - ^ i B l exp(jiZi)dzi/dai (4a) 
Expressing dz/da as (dz/dp)(dp/da) and using eqs. (3b) and (3c), one gets 

Qi = 2(N APi /A) 1 / 3 [A(C i + B t exp( 74 z 2 ) - B tll { Pt / p*) 1 ' 3 
(1 - oti ) exp(7,Zi)] 

Finally, the bulk modulus of the i-th phase of a material having the atomic mass A and density p is given by 

B = 2p l (N A e 2 /9A)(NAP l /A) 1 / 3 [A(C l + B i exp( 7i ^)- 

5,7,(1 ~ ar V8 ) _1 U/Pi) 1/8 «p(7i^)] 
The pressure derivative of B can be calculated from eqs. (1) and (2). It thus follows that 

B' = dB/dP = (d/dP)(p(dP/dp)) = d/dP[pN A e 2 Q/9A] = 
= 1 + P (N A e 2 /9A)dQ/dP = 1 + ( P /Q)dQ/dp 

and 

B' = 1 + 8C z N\ /3 /(24aN\ /3 + < 2 >)+ < 3 > +... (8) 
where C2> and < 3 > denote the number of terms omitted due to space limitations. 

Conclusions 

In this note we have briefly presented a theoretical algorithm for the determination of the isohermal 
bullk modulus and its first pressure derivative within a particular semi-classical theory of dense matter. This 
is a distinct advantage over the existing literature in the high-pressure field, where B and B' are usually 
derived by fitting the experimental p-V data to some previously chosen form of the EOS. 

The results obtained in this note give the possibility of refining the models of the internal structure of 
cold dense astrophysical objects (such as planets satellites or asteroids) as well as theoretically following the 
changes of compressibility of laboratory specimen undergoing phase transitions under high pressure. Details 
will be published elsewhere. 
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